Enhancing Constraint Propagation with 
Composition Operators * 

Laurent Granvilliers and Eric Monfroy 

Computer Science Research Institute (IRIN) 

University of Nantes 

44322 Nantes cedex 3, France 

o 
o 

(N 

Abstract. Constraint propagation is a general algorithmic approach for 

3 , pruning the search space of a CSP. In a uniform way, K. R. Apt ||l| has de- 

' ' fined a computation as an iteration of reduction functions over a domain. 

Cm I In |Q] he has also demonstrated the need for integrating static properties 

of reduction functions (commutativity and semi-commutativity) to de- 
sign specialized algorithms such as ACS and DAC. We introduce here a set 
■^C , of operators for modeling compositions of reduction functions. Two of the 

major goals are to tackle parallel computations, and dynamic behaviours 
(such as slow convergence). 
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^ . 1 Introduction 

> 
(N . 

f^ , A Constraint Satisfaction Problem (CSP) is defined by a set of domains, a set of 

^^ ' variables, and a set of constraints. Solving a CSP consists in finding assignments 

C^^ , of variables (i.e., values from their domains) that satisfy constraints. Since this 

problem is NP-hard, preprocessing techniques have been implemented to prune 
domains (search space) before backtracking, e.g.. filtering algorithms based upon 
consistency properties of subsets of constraints |10| . Constraint propagation is a 
Q ', generic term for these techniques. 

Recently, K. R. Apt Q has proposed a unified framework for constraint 
propagation. The solving process is defined as a chaotic iteration Q, i.e., an 
J^ ' iteration of reduction functions over domains. Under well-chosen properties of 

JH , domains (e.g., partial ordering, well-foundcdness) and functions (e.g., mono- 

tonicity, contractance) , iteration-based algorithms are shown to be finite and 
confluent. Further refinements have been devised in g to tackle strategies based 
on additional properties of functions (such as AC3 with idempotence and commu- 
tativity). Hence, specializations of component functions of the generic iteration 
algorithm have been implemented to tune the order of applications of functions. 
In this paper we introduce new properties of functions (i.e., strongness, in- 
dependency, and redundancy) to tackle parallel computations |q| , and dynamic 
behaviours such as slow convergence B. Furthermore, the dynamic nature of 

* This work has been supported by the 1ST project COCONUT from the European 
Community. An extended version of this paper will be submitted to the CP'2001 
conference. E-mails: {granvilliers, monfroy }@irin. univ-nantes.fr 



such strategies led us to define tlie notion of composition operator. Basically, a 
composition operator is a local strategy for applying functions: it implements 
combinators for sequential or parallel computation, and for computation of clo- 
sures. An iteration (preserving the semantics — termination, confluence — ) is 
then defined as a sequence of application of composition operators. This ap- 
proach provides flexibility (deflnition of components for creating composition 
operators) and generality (single iteration algorithm). 

Using composition operators, we then present several strategies for mod- 
eling heuristics and properties of well-known constraint solvers: priorities of 
constraints in the finite domain solver Choco Ig], slow convergence arising in 
interval narrowing ||9(], efficient heuristics for interval narrowing g as is done 
in Numerica ||l2| , and finally, parallel computations |^. Essentially, each strat- 
egy dynamically creates and applies composition operators until a fixed-point is 
computed. 

The outline of this paper is the following. A basic constraint propagation 
framework is introduced in Section |2|. Our new framework based on composition 
operators is described in Section pT In Section 0, some existing strategies are 
shown to fit in our model. 



2 Constraint Propagation 

Our aim is not to be as general as in pi . We consider here contracting and mono- 
tonic functions on a finite semilattice ordered by set inclusion. Such a domain 
is used to address the decoupling of reductions arising in parallel computations 
(see Section pi). Note that the results given here hold under these assumptions. 
The computation domain is a finite semilattice (2?, C,n), i.e., a partially 
ordered set {T>,C) in which every nonempty finite subset has a greatest lower 
bound (an inf element). The ordering C corresponds to set inclusion. The meet 
operation is the set intersection n, and the bottom element is the empty set. 
Note that in most constraint solvers, the computation domain is a finite lattice, 
i.e., finite domains such as Booleans, Integers, fioating-point intervals. 

Definition 1 (Reduction function). Consider a function f on V. 

— f is contracting ifix £ V, f{x) C x; 

— f is monotonic ifyx,y G P, a; C y => fix) ^ /(j/)- 

A reduction function on V is a contracting and monotonic function on V. 

In the following we consider a finite set of reduction functions F — {fi, . . . , fk} 
on P. 

Definition 2 (Iteration). Given an element d ^ V, an iteration of F over d 
is an infinite sequence of values do,di, . . . defined inductively by: 

do := d 

di ■= l3,{di-i) i > 1 



where ji is an element of [1, ... ,k]. A .sequence da D di ^ ■ ■ ■ of elements from 
T) stabilizes at e if for some j ^ we have di = e for i'^ j . 

We have the foUowing stabiHzation lemma. 

Lemma 1. Suppose that an iteration of F over d stabilizes at a common fixed- 
point e of the functions from F . Then, e is the greatest common fixed-point of 
the functions from F that is included in d, i.e., (n*L]^/i) t oj{d). 

Proof. See /@/. It follows from the monotonicity of the reduction functions. 

The generic iteration algorithm (Gl) defined by K. R. Apt in ||^ is given in 
Table |l|. The correctness of GI is stated by Theorem ^ 



Table 1. Generic Iteration Algorithm based on Reduction Functions. 



GI [F: set of reduction functions onV ; d: element oiV): V 
begin 

G := F 

while G 7^ do 
choose g £ G 
G := G-{g} 
G := GU update (G, g, d) 
d ■.= g{d) 
od 

return d 
end 

where for all G, g, d the set of functions update (G, g, d) is such that 

A. {/ e F - G 1 /(d) = d A /(7(d) ^ g{d)} C update (G, g, d) 

B. g[d) = d implies that update (G, g,d) = 

C. gg{d) 7^ g{d) implies that g € update {G,g,d) 



Theorem 1. Every execution of GI terminates and computes in d the greatest 
common fixed-point of the functions from F. 

Proof. See ^. It follows from the monotonicity and contractance of the reduction 
functions, and the well-foundedness of the ordering of the semilattice. 



3 Composition Operators 

We restrict our attention to simple domains. The generalization of our results to 
compound domains, i.e., considering fc-ary functions on (2?i, . . . , 2?fe), is straight- 
forward. 



3.1 Definitions 

Consider a finite set of functions F — {/i, . . . , fk} on 2?. We introduce a set of 
composition operations on F as follows: 



Sequence 

Closure 

Decoupling 



F° denotes the function 
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Note that F is supposed to be ordered since the sequence operation is not com- 
mutative. This assumption is no longer necessary for the closure and decoupling 
operations since the computation of a fixed-point is a declarative property, and 
the intersection operation is commutative. 

There are several motivations for introducing such operations in a generic 
propagation framework: 

— The sequence operation fixes the order of application of the reduction func- 
tions. It can be used for computing directional arc consistency based on the 
semi-commutativity property, for modeling priorities of solvers, and for im- 
plementing heuristics or knowledge of solvers about their relative efficiencies. 

— The closure operation allows us to make a closure from a non idempotent 
function, and to describe multi-level algorithms that compute fixed-points 
of different solvers. 

— The decoupling operation is essentially used to model parallel computations, 
enforcing different functions over the same domain, and then computing the 
intersection using a fold reduction step. 

The notion of composition operator models a function (a complex solver) 
built from composition operations. 

Definition 3 (Composition operator). Let F = {/i, . . . , fk} be a finite set of 
reduction functions. A composition operator on F is a function V —> V defined 
by induction as follows: 

Atomic: fi is a composition operator for i — 1, . . . ,k. 

Given a finite set (p of composition operators on F, 



Sequence 

Closure 

Decoupling 



<?° is a composition operator; 
<?" is a composition operator; 
'P'^ is a composition operator. 



The generator Gen(0) of a composition operator (p on F is the subset of functions 
from F that are "involved" in (j). It is defined inductively by: 

— Gen(0) = {(j)} if (j) is an atomic operator from, F; 

- Gen((^) =Gen((/ii)U---UGen((/>fc) i/ = {0i, . . . , 04* /or • G {o,LU,n}. 

Lemma states that a composition operator is also a reduction function. 



Lemma 2. Consider a finite set of reduction functions F ^ {fi, . . . , fk}- Then, 
every composition operator on F is 

(i) contracting; 
(a) monotonic. 

Proof. The proof is done by induction. Every atomic operator is contracting 
and monotonic by definition of a reduction function. Now consider a set <l> = 
{01, . . . ,4>k} of contracting and monotonic composition operators. 

(i) By hypothesis the composition operators from <P are contracting. Then it is 

immediate to prove the contractance of <1>° , <1>^ , and (p'^ . 
(a) Given x,y (z V suppose that a; C y. 

- Since every (pi is supposed to be monotonic, then, we have 4>k{x) C (j)k(y), 
and then, (/jk-itpkix) C 0j,_i(/)fc(y), and so on. It follows that<P° is mono- 
tonic. 

- To prove the monotonicity of 'P'^ we consider the function ip : n^^it/)^. 
Then, we prove that ip is monotonic (third item). It follows that ip f 
u}(x) C (^ f a;(y). Then it is immediate to prove (by a double inclusion) 
that the set of fixed-points of p) coincides with the set of common fixed- 
points of the functions from <I>. As a consequence, we have 1^9 f w = ^", 
that completes the proof. 

- For i — 1, . . . ,k, we have (f>i{x) C (j)i{y) by monotonicity of 4>i. It follows 
that n*L]^(/),;(x) C CXl^i4>i{y) , that ends the proof. 

Lemma states that a fixed-point of a composition operator is a common 
fixed-point of the functions from its generator. The key idea is that the ap- 
plication of a composition operator imphes that each reduction function in its 
generator is apphed at least once. 

Lemma 3. Consider a finite set of reduction functions F and a composition 
operator (j) on F . Then, e is a fixed-point of (f> if and only if e is a common 
fixed-point of the functions from Gen((/)). 

Proof. We prove by induction the equivalence <j){e) = e 4^ V/ G Gen((/)) /(e) = e. 
It obviously holds for an atomic operator (j). Now consider a set of composition 
operators <I> — {(j)i, . . . , 0/c} on F and assume that the equivalence holds for each 
4>i, . . . ,4'k. If we have (/) = $* for -k £ {o, uj, n}, then it follows: 

(j){e) = e <^ Vi G {1, . . . , /c} (f>i{e) — e (immediate result) 

<^ Vi G {1, . . . ,k} \ff G Gen(0i) /(e) = e (induction hypothesis) 

Considering Gen(0) = UiGen{(/)i) completes the proof. 

The notion of iteration is slightly extended to deal with composition operators 
instead of reduction functions. 



Definition 4 (Iteration). Consider a finite set of reduction functions F, and 
a finite set (l> = {(j)i, . . . ,(f)k} of composition operators on F. Given an element 
d GV, an iteration of<P over d is an infinite sequence of values do, di, . . . defined 
inductively by: 

do := d 

d, := 4>j^{di-i) i > 1 

where ji is an element of [1, ... ,k]. 

Lemma n^ (stabilization lemm.a from K. R. Apt) remains valid, that follows 
from the monotonicity of the composition operators (Lemma 0). Moreover, we 
have the following, essential result. 

Lemma 4. // an iteration on <P — {(jji, . . . , cjjk} over d stabilizes at a common 
fixed-point e of the functions from <1>, and F — UiGen((j)i), then e = F'^(d). 

Proof. We first prove that e is a common fixed-point of the functions from F. 
By hypothesis e is a fixed-point of each (pi € <1>. By Lemma U it follows that 
/(e) — e for each f e Gen(0i). The proof is completed since by hypothesis the 
set of generators covers F . 

We prove now that e is the greatest common fixed-point of the functions from 
F . Consider a common fixed-point e' of the functions from F . It suffices to 
prove that e' is included in every element from the iteration, namely do, di, . . . . 
It obviously holds for i = 0. Suppose now it holds for i, i.e., e' C di, and assume 
that di^i — 4>j{di) for some j € [1,...,A:]. By monotonicity of (f)j we have 
4'j{^') ^ <i>j{di). By Lemma^ we have 4>j(e') — e' , that completes the proof. 



Table 2. Generic Iteration Algorithm based on Composition Operators. 



GIco(_F: set of reduction functions on © ; d: element oi D): T) 
begin 

G := F 

while G / do 

4> := create a composition operator on G 
G := G-Gen(0) 
G : = G U update (G, (j>, d) 
d := (j){d) 
od 

return d 
end 

where for all G, 0, d the set of functions update (G, 4>, d) is such that 

A. upA := {/ € F - G I /(d) = d A /<^(d) / (t>{d)} C update (G, 0, d) 

B. upB := (j>{d) = d implies that update (G, (ji, d) = 

C. upC := {/ G Gen(</>) | /(^(d) / (^(d)} C update (G, <f), d) 



We describe now a generic iteration algorithm GIco based on composition 
operators on a finite set of reduction functions F (see Table 0). Note that the 
set of composition operators is not fixed, since each operator is dynamically 
created from the set G of active reduction functions, and it is applied only once. 
Nevertheless, Theorems proves the correctness of GIco with respect to F. 

Theorem 2. Every execution of GIco terminates and computes in d the greatest 
common fixed-point of the functions from F. 

Proof. The proof is a direct adaptation of Apt's ||/. To prove termination it 
suffices to prove that the pair {d, 4f^G) strictly decreases in some sense at each 
iteration of the while loop, and to note that the ordering C is well-founded. 

The correctness is implied by the invariant of the while loop, i.e., every 
f G F — G is such that f{d) = d. It follows that the final domain is a common 
fixed-point of the functions from F (since G = 0). The second part of the proof 
of Lemma ensures that it is the greatest one included in the initial domain. 

The following corollary concerns the application of a closure operator in 
algorithm GIco. 

Corollary 1. Consider operator (j) that is applied in algorithm GIco. // is 
idempotent, then assumption C is reduced to upC := 0. 

Proof See Apt g/. 

3.2 On Properties of Functions 

In this section we examine the properties of reduction functions that are pre- 
served by compositions. We first define a set of properties of interest for our 
purpose. 

Definition 5. Consider two reduction functions /, g on T>. Then, for all x £ T>: 

— f is idempotent if ff{x) ~ f{x) 

— f and g commute if fg{x) — gf{x) 

— f semi- commutes with g if gf{x) C fg{x) 

— f is stronger than g if fix) C g[x) 

— f and g are independent if fg{x) ~ gf{x) = f{x) r\g{x) 

— f and g are redundant if f{x) = g{x) 

— f and g are weakly redundant if f 1 uj{x) — g 1 uj{x) 

The property of idempotence of a composition operator allows us to modify 
assumption C on the update function of the GIco algorithm (see Corollary |l|). For 
this purpose some kinds of composition operators are shown to be idempotent 
in the following proposition. 

Proposition 1. Consider a set of composition operators (p — {0i, . . . , 0^}. 
Then, (j) is idempotent if: 

i) 4> is a closure operator <P'^ . 



a) (j) is a sequence operator <P° , (pi semi-commutes with (f)j for i > j , and each 

<j)i is idempotent. 
Hi) (j) is a sequence operator (p° , for each i there exists j < i such that (j)j is 

stronger than (j)i, and (pi is idempotent. 
iv) (j) is a decoupling operator (p'^ , each (pi is idempotent, and the (pi are pairwise 

independent. 

Proof. It suffices to prove that in all these cases, (p is idempotent. The proof 
is then completed by Lemma H, i.e., each element computed by (p is a common 
fixed-point of the functions from its generator. 

i) The proof is obvious, 
ii) See Apt J^. 
Hi) Since the relation of strongness is transitive, then (pi is stronger than 4>i for 

i G [2, . . . ,k]. Now it suffices to prove that (pj(pi . . .(pk{x) = (pi . . . (pkix). We 

have 

<t>i{y) = <Pi<Pi{y) ^ <t>]<t>i{y) ^ (pi{y) 

since (pi is idempotent, (pi is stronger than (pj, and (pj is contracting. This 
ends the proof if we set y = (p2 ■ . . (pk{x). 
iv) Let us prove that function n^^j^t/)^ is idempotent. We prove by induction on j 
that C\'[^^(pi is idempotent for j = I, . . . ,k and that it is independent on (pj+i 
for j — 1, ... ,k — 1. It holds for j ~ 1 since by hypothesis, (pi is idempotent, 
and (pi and (p2 are independent. Now fix 1 < j < k, and consider that 
Lp : C\l^-^(pi is idempotent, and that ip and (pj+i are independent. We prove 
that function ip : ipr\(pj^i is idempotent and independent on (pj+2- Given an 
element x G T>, we have: 

i^ii^ix)) = (p{(p{x) n (pj+i{x)) n (pj+i{ip{x) n (pj+i{x)) 

= (pip{(l)j+i{x)) n (pj+i(pj+i{(p{x)) independence ip,(pj+i 

= (p(pj+i{x) n (pj+iip{x) idempotence tp, (pj+i 

= {f(x) n (pj+i{x)) n {(pj+i{x) n '^{x)) independence ip, (pj+i 

= ^j{x) commutativity of n 

Then ip is idempotent. Now we prove that ip and (pj+2 o.i'e independent, i.e., 

{ip n (pj+i){(l)j+2{x)) = (t>j+2{i'f n (pj+i){x)) == ((y3 n (j)j + l){x) n (pj+2{x) 

It suffices to prove that each term of this formula is equivalent to H^^Z^j^ (x). 
It obviously holds for the last term. For the first term, we have: 

{ip n (pj+i){(j)j+2{x)) = Lp(j)j+2{x) n (pj + i(pj+2{x) 

= (pi(pj+2ix) n • • • n (pj(pj+2ix) n 0^+10^+2(2;) 

The independence of each pair {(pi, (pj) ends the proof. For the second term 
we have: 

4>J+2{{^ n (pj + i){x)) = (pj+2ip{x) n 0j+20j + l(x) 



It suffices to remark that (j)j+2 cind '/>j+i are independent by hypothesis, and 
then to prove that (t>j+2 o,nd ip commute, i.e., (j)j+2V{x) = Lp(t)j+2{x). This 
result is easily proved by induction since 

(f>j+2(p{x) = 0J+201 • ■ ■ (t)jix) 

by independence of each pair {(pi, cpj), and then 

0J+201 • • ■ 4>j{x) = 4>l4>j+24'2 ■ ■ ■ 4>j{x) 

by independence of (J3j+2 o,nd (j)i, and so on. 

In the following proposition, we identify cases where the computation of a 
closure of a set of composition operators can be improved, according to properties 
of independence and redundancy of operators. 

Proposition 2. Consider a set of composition operator <l> and a composition 
operator ip. Then, the following properties hold: 

1. If for all (f) ^ <P, ip and cj) are independent, then {<P U {p})" ~ 'T'^ f^ p ] uj. 

2. If there exists (j) G <P such that p and <j) are weakly redundant, then (<P U 

Proof. 1. It suffices to show that operator 9 : (P^ C\ p ] lo is idempotent, i.e., 
99{x) — 9{x) for all x G T). The proof then follows by Lemma Q. Given a 
particular element x, consider that 9{x) = 4>ii ■ ■ ■ (pik (a;)n(/3 . . . p{x), and that 
99 (x) — (f)j-^ . . . (f>j^ 9(x)np . . . p9(x) . A simple induction, using the hypothesis 
of independence, then allows us to rewrite 99{x) as 9{x) fl 9{x). 
2. The proof is obvious since by definition, <I>^ is a fixed-point of 4>. Hence, it 
is also a fixed-point of p. 

3.3 Prom Decomposition to Composition 

A glass-box solver mainly combines elementary solving components. The relation 
between solvers and components is expressed by Definition pi 

Definition 6 (Decomposition of a function). Consider a reduction function 
f on T>. A finite set of reduction functions F on T) is a decomposition of f if 
fj^F > 1 and for every x £ V, we have F^{x) — f 1 uj(x). 

For instance consider a filtering algorithm enforcing a local consistency technique 
over a CSP. The generic computation is an iteration such that the consistency 
of a set of constraints, a constraint, or a constraint projection is verified at 
each step. Note that there is a correspondence between the decomposition of a 
data (e.g., a CSP) in a set of elementary components (e.g., constraints), and 
the decomposition of a function (e.g., a solver associated with a CSP) in a 
conjunction of elementary functions (e.g., an elementary solver associated with 
constraints) . 

Considering different levels of granularity in a constraint solving process must 
not infiuence the semantics of computation. To this end we have the following 
result. 



Proposition 3. Consider a finite set of reduction functions F = {/o, /i, . . . , fk} 
on v. Let G be a decomposition of /q. Given H = G U {fi, . . . , fk}, we have 
H'^{x) = F'^{x) for every x eV. 

Proof. The proof is very similar to the one of Lemma 0. 

The decomposition of reduction functions leads to the notion of decompo- 
sition relation. It can be useful to compare levels of granularity of iteration 
algorithms. 

Definition 7 (Decomposition relation). Consider F,G two finite sets of re- 
duction functions on V. G is a decomposition of F if ^G > i^F and G'^{d) = 
F'^{d) for every de V. 

Note that the decomposition relation is a strict partial order. 

In practice the aim is to fastly design efficient algorithms. We believe that 
these requirements can be achieved in a glass-box and generic (e.g., object ori- 
ented) programming approach. 

— Developing powerful propagation techniques can be tackled by the GIco algo- 
rithm taking as input a well-chosen decomposition of solvers. The generation 
of composition operators during the iteration allows us to efficiently com- 
bine elementary solvers using the best strategy with respect to the knowledge 
or properties of reduction functions. Several efficient strategies, dynamic in 
essence, are described in Section Q. 

— Efficient prototyping is achieved by a generic programming approach. In fact 
in the case a generic (constraint propagation) algorithm is re-used, only a 
(generic) piece of code for each kind of elementary solver has to be imple- 
mented, e.g., for a reduction function enforcing a local consistency property. 
In practice a main challenge is then to automatically generate the set of 
reduction functions of the decomposition given a generic implementation of 
a function and a CSP. This issue is discussed in the conclusion but a precise 
description is out of the scope of this paper. 



4 Specialized Algorithms 

The constraint propagation framework described in this paper allows us to tackle 
existing, efficient algorithms based on strategies and heuristics. 

4.1 Scheduling of Constraints with Priorities 

Choco is a constraint programming system for finite domains that has been 
developed by Laburthe et al. g]. The core algorithm is constraint propagation 
whose main feature is to process constraints according to the complexity in 
time of their associated solving algorithms. Thus, primitive constraints are first 
processed, and then, global constraints with linear, quadratic complexity, and so 
on until a fixed-point is reached. In other words the computation is a sequence 



of application of closures, each closure being connected to the previous one by 
means of propagation events (modification of domains). 

Implementing the propagation engine of Choco using our generic algorithm 
GIco can be done by considering priorities of reduction functions. The creation 
of composition operators can be implemented in two ways as follows: 

1. 4': = {g G G I priority{g) — a}'^ s.t. a = niin( {priority (5) | g G G}) 

2. 0:=(G^U---UG^)° s.t. Gi,...,Gp is a partition of G 



Vi e [1, . . . ,p],Vg e Gj,priority(5) ^ 
ap < ap^i < • • • < ai 



a,; 



The first implementation describes the computation of the closure of the set of 
active reduction functions with the greatest priority (i.e., the computationally 
less expensive functions). It is a model of Choco in the sense that only functions 
with greatest priority are applied at a time. Note that the composition operator 
is a closure (and thus, idempotent), and consequently Corollary |l] applies for the 
update function (i.e., upC= 0). 

The second implementation is a sequence of closures, each closure processing 
the set of functions of a given priority. The main difference with respect to the 
first method is that no propagation step is performed between the application of 
two closures. This approach tends to minimize costs for updating propagation 
structures. 



4.2 Sequentiality in Interval Constraints 

The solver presented in g implements constraint propagation for interval do- 
mains, where reduction functions enforce box consistency for numeric constraints 
over the reals. It extends two existing solvers, namely Numerica Jl^ by Van Hen- 
tenryck et al., and the algorithm BC4 0] by Benhamou et al. 

The solving process combines three kinds of reduction functions: 

1. Function fc^i computes box consistency for the domain of variable Xi with 
respect to constraint c, i.e., /c,i is a projection function; 

2. Function gc implements constraint inversion for all variables of constraint c; 

3. Function h computes a linear relaxation (by means of a first order Taylor ap- 
proximation) for a constraint system, which is then processed by the interval 
Gauss-Seidel method. 

The best strategy is based upon properties and heuristics: function fd is stronger 
than the projection of gc on Xi, but they are redundant if Xi occurs only once in 
c; gc is computationaly less expensive than functions /c,i for all i; h is in general 
more precise for tight domains, while the fd and gc are more efficient for large 
domains. Using GIco this strategy is efficiently implemented as follows: 

ci,:={{G-{h}ru{hrr 

Note that h depends on all variables, that implies that h belongs to G. 



Furthermore, the decomposition process that generates the set of reduction 
functions is tuned with respect to the redundancy property. Hence, each re- 
duction function fc^i such that Xi occurs only once in c is removed (since gc is 
as precise as fc^i for Xi). Proposition guarantees that the output domain is 
consistent with respect to all constraints from the initial system. 

4.3 Acceleration of Interval Narrowing 

Interval narrowing, i.e., constraint propagation with interval domains, is ineffi- 
cient if slow convergence happens. A slow convergence corresponds to a cycle of 
reduction functions fi . . . fj . . . fi such that each application only deletes a small 
part of a domain. When constraints are nonlinear constraints over the reals, this 
problem frequently arises, due to, e.g., singularities or points of tangence. 

Lhomme et al. M have devised an efficient strategy based on cycle detection 
and simplification. The aim is to locally select and apply the best reduction 
functions while delaying some active functions supposed to slow the computation. 
Given a cycle, i.e., a set of functions G' from the propagation set G, the solving 
process using GIco can be described as follows: 

1. G" is rewritten as <? U <?' where <P' contains all functions g from G such that 
for all (f) £ G — {g}, g and (/) are independent; 

2. ^ is rewritten as <?i U <?2, ^i being composed of the best functions from <P. 
More precisely 'Pi contains a function / per variable whose (current) domain 
can be modified by a function from ^; the selected function / is the one that 
computes the largest reduction; 

3. Doing so, the composition operator applied in GIco can be defined as: 

applying first the best functions, and then the ones that have been delayed 
because of the independency property. Note that each function from the set 
(G — G') U ^2 has to be added in G by the update function. 

4.4 Parallel Constraint Propagation 

Parallel processing of numerical problems via interval constraints has been pro- 
posed as a general framework for high-performance numerical computation in H . 
Parallel constraint propagation mU operationally consists in distributing reduc- 
tion functions among processors, performing local computations, and then ac- 
cumulating and intersecting new domains on some processors. 

The decoupling composition operator can be used to implement parallel con- 
straint propagation. A basic strategy is to create a partition of the propagation 
structure G = Gi U • • • U Gk, k being dependent on the number of processors, 
and to consider operator (/) to be applied in the algorithm GIco: 

0: = G°n...nG^ 



Moreover if one wants to perform more local computations before synchronisa- 
tion and communication, a closure can be computed on each processor as follows: 

Nevertheless, it has been observed that the classical notion of parallel speed- 
up is not a correct measure of success for such algorithms. This is due to a 
parallel decoupling phenomenon: convergence may be faster when two interval 
contractions are applied in sequence than in parallel. As a consequence a par- 
allel version of Lhomme's strategy, described in the previous section, has been 
proposed in Q . Essentially, parallelism is only used to select the best functions, 
i.e., to create the decomposition ^i U <?2 of the previous section. 

Note that the decoupling phenomenon can be controlled according to the in- 
dependence property of reduction functions. Proposition S ensures that a closure 
can be computed in parallel, provided that each subset of dependent functions 
is located on one processor. This property is not achievable in general for the 
whole CSP. Nevertheless the number of links between processors (corresponding 
to couples of dependent functions) can be minimized, which tends to maximize 
the amount of reductions of domains. In that case it may also be efficient to 
duplicate some functions on several processors in order to break some links. Fur- 
ther work will investigate such strategies, which have to be dynamic to guarantee 
load balancing. 

5 Conclusion and Perspectives 

A set of composition operations of reduction functions is introduced to design 
dynamic constraint propagation strategies. K. R. Apt's iteration model is slightly 
modified while preserving the semantics. Finally, several well-known strategies 
(using priorities of constraints, heuristics on the order of application of functions, 
and parallelism) are modeled using a single iteration algorithm. 

A generic implementation of constraint propagation, integrating composition 
operators, has been designed. However it is out of scope of this article and it will 
be the topic of a second article. 

The set of composition operators is (intentionally) reduced to sequence, clo- 
sure, and decoupling operators. One may desire additional operators to model 
sequences of fixed length, quasi closures with a notion of precision, or conditional 
strategies with respect to dynamic criteria. We believe that their integration in 
our framework is feasible. 
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